Origin of the Efficient Polaron Pair Dissociation in Polymer-Fullerene Blends 
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The separation of photogenerated polaron pairs in organic bulk heterojunction solar cells is the intermediate 
but crucial step between exciton dissociation and charge transport to the electrodes. In state-of-the-art de- 
vices, above 80% of all polaron pairs are separated at fields of below lO^V/m. In contrast, considering just the 
Coulomb binding of the polaron pair, electric fields above lO^V/m would be needed to reach similar yields. In 
order to resolve this discrepancy, we performed kinetic Monte Carlo simulations of polaron pair dissociation 
in donor-acceptor blends, considering delocalised charge carriers along conjugated polymer chain segments. 
We show that the resulting fast local charge carrier transport can indeed explain the high experimental quantum 
yields in polymer solar cells. 
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Polymer-fullerene solar cells have seen steady improve- 
ments in terms of performance in the course of the last years. 
Bulk heterojunction devices showing 5-6% power conversion 
efficiency have been published last year,' which is only pos- 
sible by achieving internal photon conversion efficiencies of 
beyond 80%. The light is mainly absorbed in the conjugated 
polymer, generating singlet excitons, which diffuse towardj^ 
the interface to the electron accepting fuUerene derivative.^ 
The ultrafast electron transfer to the acceptor molecule is very 
efficient.^ Electron and hole are now spatially separated but 
still Coulomb bound. Akeady a few years ago, the importance 
of optimizing the morphology of the donor-acceptor phases in 
these bulk heterojunction devices for an efficient polaron pair 
separation was shown empirically. A major breakthrough to 
achieve this was the selection of a suitable solvent.^ However, 
even with an optimized phase separation, it was the common 
understanding that the dissociation was mainly driven by the 
external electric field. In this paper, the origin of the efficient 
polaron pair separation in bulk heterojunction solar cells is re- 
considered. 

The dissociation of polaron pairs is generally expected only 
when electric fields exceeding lO^V/m are applied.l^In a bulk 
heterojunction solar cell, however, the internal field in the rel- 
evant fourth quadrant of the cuiTent-voltage characteristics 
reaches its maximum value of often less than lO^V/m un- 
der short circuit conditions. Despite the electric field being 
more than one order of magnitude smaller, internal quantum 
efficiency of 80% under short circuit are determined experi- 
mentally. This large discrepancy is particularly surprising as 
the polaron pair lifetime was recently determined to be on the 
order of only lOns,^ seemingly making the polaron pair disso- 
ciation even more unfavorable and challenging. As the classi- 
cal models such as Onsager^ and Braun** cannot predict these 
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high yields with Ufetimes so shoit, effective lifetimes in the 
microsecond range have typically been used in modelling and 
macroscopic simulations.!^ 

In order to find the additional driving force assisting the 
geminate pair separation, Arkhipov et al. proposed the ana- 
lytic dark dipole model, based on the experimental finding of 
vacuum level shifts at the heterojunction of oligomer/fuUerene 
bilayers."* We point out that optical investigations on spin 
coated polymer films show that the polymer chains were pref- 
erentially oriented in parallel to the substrate, with dipole mo- 
ments laying almost parallel to the chains. This effect has 
been shown to be proportional to the molecular weight,'^ in- 
dicating that a remaining net dipole moment is expected to be 
present in many conjugated polymers. Another reporP ap- 
plies a Monte Carlo simulation to show that a large initial 
separation of the constituents of polaron pairs yields higher 
dissociation yields, although an additional space charge re- 
gion had to be assumed to explain the experimental photocur- 
rents of bilayer solar cells. Hot polaron pairs being formed 
after exciton dissociation, with excess kinetic energy avail- 
able for the dissociation process, were proposed as model pic- 
ture. A further approach to improve the dissociation is the 
high local mobility within ordered domains in the otherwise 
disordered donor-acceptor blends, as studied for instance by 
Groves et al!^ However, while the high experimental yield is 
almost reached, this is due to using large donor and acceptor 
grains, and polaron pair recombination rates orders of magni- 
tude lower than determined experimentally.'^ 

In this paper, we present kinetic Monte Carlo simulations of 
hopping transport of Coulomb bound polarons in an energeti- 
cally disordered donor-acceptor blend system. We investigate 
the geminate pair dissociation, and how it is influenced by a 
charge carrier delocalization due to extended polymer chains. 
For the first time, it is shown that the very efficient polaron 
pair separation responsible for the high performance of state- 
of-the-art organic solar cells can be explained by conjugation 
lengths in the range of 7-10 units. 
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The simulation space was a simplified bulk heteroj unction 
system, a simple cubic lattice of 100 units length sandwiched 
between two electrodes for carrier extraction. The lattice, with 
a constant spacing of Inm, had an area of 30 x 30 units with 
periodic boundary conditions. A model for extended, hole 
transporting polymer chains — similar to the slithering snake 
approach in Ref.'"*^ — was implemented. The charge carriers 
within a single segment (the effective conjugation length CL 
of the polymer chain), given by a fixed number of monomer 
units, were assumed to be delocalized. For instance, for a 
polymer chain consisting of 4 monomer units, CL — 4, delo- 
calization was implemented by placing partial charges (here, 
quarter charges) on each of these four units. According to the 
Inm distance between neighboring sites, this segment is 4nm 
long, although mostly not straightened. To leave this polymer 
chain segment, the whole charge can hop 'as one' to the next 
segment, where it spreads again. As the charge transport on 
a segment (intra-chain transport) is much faster than the hop- 
ping process between chains (inter-chain hopping), the former 
is assumed to be instantaneous. Charge transport between seg- 
ments was implemented using the Miller-Abrahams hopping 
rate, 

v,=Voexp(-2y.,)|^"P("^) '^^'^^^ (1) 

where Vq = lO'^s^' is the attempt-to-escape frequency and 
Y==5- lOV-i the inverse localization radius. r,y denotes the 
hopping distance, and AEij the energy difference of site i and 
site j, also accounting for the external electric field and the 
Coulomb force between all charge carriers and mirror charges. 
Alternatively to hopping rate Eqn.[T[ also a polaronic hopping 
rate could have been used,"^ which might change minor de- 
tails of the simulation. Half of all sites was defined as elec- 
tron accepting and transporting molecules, such as fuUerene 
derivatives, each occupying just a single site. Energetic disor- 
der of the density of states was assumed, with a width of the 
Gaussian distribution of 75meV for the polymer donor and 
60meV for the fullerene acceptor, according to our measure- 
ments on poly(3-hexyl thiophene) (P3HT) and [6,6]-phenyl- 
Cgi butyric acid methyl ester (PCBM).ESI We considered ten 
polaron pairs simultaneously, each consisting of one positive 
polaron on a donor site and one negative polaron on an accep- 
tor site, yielding a carrier concentration of about lO^^cm^^. 

At least 200 runs were performed for each parameter set. 
We consider a polaron pair separated if at least one con- 
stituent is extracted at an electrode, thus escaping geminate 
and nongeminate recombination. In addition to considering 
CL, we replaced the polaron pair lifetime T/^ by an effective 
'Ceff > '^f- Thus, we can account for a reduced recombination 
probability potentially induced by additional effects such as 
interface dipoles, which are of theoreticaff^l and experimental 
relevance. ^ ' 

The influence of the effective lifetime on the polaron pair 
separation yields at a constant donor conjugation length of one 
is presented in Fig.[TJa). A lower recombination rate leads to 
a higher separation yield. Meanwhile, even at effective life- 
times as long as 1ms, the high experimental dissociation effi- 
ciency at low electric fields cannot be explained. The variation 
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FIG. 1: (Color Online) Polaron pair dissociation yield in 1:1 donor- 
acceptor blends at 300K. The yellow rectangle denotes the typical 
internal field under short-circuit conditions of organic solar cells, (a) 
Variation of the effective lifetime tgff at a constant donor CL of 1. 
The efficient experimental dissociation at low electric fields cannot 
be reached (dielectric constant e = 3.0). (b) Variation of the conjuga- 
tion length CL, with teff fixed at ljUS (dielectric constant e = 3.0). (c) 
Simultaneous variation of Teff and CL, with the dielectric constant 
chosen as 3.7. For long polymer chain segments, the high experi- 
mentally observed yields are replicated. 

of CL, and its influence on the geminate pair dissociation, is 
shown in Fig. [TJb) for a fixed effective lifetime of 1/js. As 
expected, longer chains lead to an improved separation yield, 
and show a more positive impact on the yield at low fields as 
compared to the reduced recombination rate. A longer conju- 
gation length of the donor leads to an increasing local charge 
carrier mobility and a larger initial electron-hole distance, as 
the positive charge is delocalised within the conjugated seg- 
ment. Therefore, the hole can escape the mutual Coulomb 
attraction more easily, and also needs less hops during its es- 
cape. 

The influence of a simultaneous variation of effective life- 
time and CL is shown in Fig.[TJc). We point out that the typ- 
ical conjugation length in polymers used in organic bulk het- 
erojunction solar cells, such as regioregular poly(3-hexyl thio- 
phene), is around 7-10 monomer units. For this regime, the 
highly efficient polaron pair dissociation process in organic 



3 



solar cells — as expected from experiment — is for the first time 
also achieved in Monte Carlo simulations when assuming ef- 
fective lifetimes between 100ns and 10/js. The prerequisite for 
this high efficiency is the physically justified delocalisation of 
positive polarons within conjugated polymer chain segments. 

Veldman et alP experimentally determined the polaron pair 
lifetime in a polyfluorene-fullerene blend to about 4ns. Val- 
ues for the more common poly(3-hexyl thiophene)-fullerene 
mixture have to our knowledge never been determined. The 
remaining discrepancy between the effective lifetime of 100ns 
or more to the experimental lifetime Zf can be due to other 
mechanisms reducing the polaron pair recombination yield. 
We define Tgff = '^f/Prr, where pr,- < 1 denotes the reduced 
recombination probability. For Teff — IfJS and if = 4ns, 
Prr = 1 /250. Considering the Miller- Abrahams hopping rate, 
Eqn. [T] this could either correspond to an energetic barrier 
for recombination of AE — A;rin(250) = 138meV at room 
temperature — e.g., due to dark interface dipoles similar to the 
considerations of Arkhipov"* — or to a reduction of the tun- 
neling term for recombination. The latter could be explained 
by a weakened wavefunction overlap between the donor and 
acceptor molecule; a reduction of the recombination proba- 
blility by factor 250 is achieved by increasing the inverse lo- 
calisation radius by factor ln(250) = 5.5. Concerning the po- 
tential influence of dipoles at the donor-acceptor interface on 
the geminate pair separation, we found it to be not significant 
(not shown). This is due to the system being a bulk hetero- 
junction, with surface normals pointing in all directions, thus 
cancelling each other on longer ranges beyond one nm. The 
short range influence of the dipoles on sub-nm scale could 
only be considered by an increased effective lifetime, which 
shows a clear impact on the dissociation efficiency. Two 
other potential means of increasing the polaron pair dissoci- 
ation is the high interchain mobility due to lamellae as re- 
ported for poly(3-hexyl thiophene),'** or the influence of im- 
age force effects if donor and acceptor inhibit different dielec- 
tric constants.'^ Thus, although the impact of the high donor 
conjugation length on the polaron pair separation efficiency is 
dominant by far, additional effects driving the dissociation are 
plausible and expected by us. 

In order to relate our simulations to an established model, 
we calculated the separation yield by the Onsager^ -Braun** 
(OB) theory, where the dissociation probability P{F) is given 
by 
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F is the electric field, ken = ^eff the effective polaron pair re- 
combination rate to the ground state, /j the sum of electron and 
hole charge carrier mobility. The field dependent dissociation 
rate kd{F) — /jk^IF) is given by 
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q/j/ezo is the Langevin recombination factor^ rpp 

is the Coulom- 



where y 

is the initial polaron pair radius, Ei, °^ ^/^pp 
bic binding energy of the pair, kT the thermal energy, Ji the 
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FIG. 2: (Color Online) Polaron pair dissociation yield with Monte 
Carlo simulation (solid lines) and the OB model (dashed lines) in 
1 : 1 donor-acceptor blends at 300K (dielectric constant e = 3.0). 
The contours are denoted with the yield values. The yellow rectangle 
denotes again the typical internal field under short-circuit conditions. 



Bessel function of order one, and b = F / {Sneeo{kT)^) . e is 
the elementary charge, and eeo the effective dielectric constant 
of the organic semiconductor blend. 

We considered a spatial distribution of polaron pair sepa- 
rations, as suggested by the Blom group. ^' The dissociation 
yields by Monte Carlo simulation and OB model at 300K 
are shown in Fig. J2] For the OB model, we set the /JTeff- 
product to 7 • 10^'^^A^ — for instance, 10^^m^/Vs-700ns — 
assuming CL = 2.5rpp, as the polymer chains are not straight- 
ened due to entanglement. 

Comparing the results of the Monte Carlo simulation with 
the OB model, as shown in Fig. [2] we found a qualita- 
tively good agreement for dissociation yields between 20 and 
75%. The Monte Carlo simulation reaches separation yields 
of above 90% for conjugation lengths beyond 9 units even at 
low electric fields. Here, the OB theory shows its limitations, 
as it cannot describe this behaviour The reason is that the 
theory does not account for hopping transport and energetic 
or spatial disorder. 

Recently, experimental findings of the Jans sen group in 
terms of an efficient free polaron generation — performing 
photoinduced absorption measurements on blends of fluo- 
rene copolymers with PCBM — led the authors to propose ex- 
tended, nanocrystalline fullerene clusters which could yield a 
high local charge carrier mobility, orders of magnitude beyond 
the macroscopically determined values.!^ Using these high lo- 
cal mobilities in the OB model, the authors achieved good fits 
of their experimental data, assuming cluster sizes of at least 
15nm diameter. Solid-state nuclear magnetic resonance re- 
laxometry investigations^^ performed on composites of P3HT 
and PCBM, however, lead to the conclusion that the fullerene 
derivative forms amorphous, not nanocrystalline phases.Sl 
Illustrating the ongoing debate, Campoy-Quiles et al. dis- 
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agree, pointing out that PCBM can form ordered domains.'--' 
On the other hand, the extended effective conjugation length 
found in polymers is well-knownP3Thus, we propose that the 
high local on-chain mobility of the conjugated polymer^^ due 
the extended conjugation length is — at least to a significant 
degree — ^responsible for the very efficient geminate pair disso- 
ciation mechanism in organic solar cells. We point out that the 
charge transport in the potentially nanocrystalline fuUerene 
phases as well as along semicrystalline lamellae regions found 
in some polymers such as P3HT'^ can further increase the 
effect of fast interchain transport along conjugated polymer 
chains considered in our simulations. 

In conclusion, by performing kinetic Monte Carlo simu- 
lations in donor-acceptor blend systems relevant for organic 
bulk he teroj unction solar cells, we found that the highly effi- 
cient polaron pair dissociation can be explained by delocalised 
charge carriers within conjugated segments of the polymer 
chain. The resulting local charge carrier mobility is much 
larger than the macroscopic one. Together with the reduced 
Coulomb attraction due to the accordingly increased intial po- 



laron pair radius, the high on-chain mobility is essential to ex- 
plain the high polaron pair separation yield. Assuming poly- 
mer conjugation lengths between 4 and 10 monomer units, we 
could simulate 60 to 90% dissociation yield at moderate fields 
of below lO^V/m, as found in working organic solar cells. 
These values correspond to recent experimental findings. We 
qualitatively fit our simulation results with the Onsager-Braun 
model. Our results show that this model is not well suited to 
describe polaron pair dissociation at low and high fields, but 
works well in the intermediate range. We point out that the 
polaron pair separation is only weakly field dependent in the 
working regime of organic solar cells, with a finite, signifi- 
cant separation yield at zero internal field. This has important 
consequences for the modelling and optimisation of organic 
photovoltaic devices. 
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